DISCRETIZED ROTATION HAS INFINITELY MANY 
PERIODIC ORBITS 
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Abstract. For a fixed A G (—2,2), the discretized rotation on 1? is 
defined by 

(x,y) i-> (y, -[x + Xy\). 
We prove that this dynamics has infinitely many periodic orbits. 



1. Introduction 

Space discretization of dynamical systems attracted considerable interests 
among researchers [91 [27J 12H 13 [6] . One motivation is to understand the 
distance between original dynamics and its computer simulation through 
discretized model. In this paper, we are interested in a discretized planar 
rotation. It is very simple but we know surprisingly little on this discretized 
system. We start with a conjecture studied by many authors, for e.g., in 
[20\ [25l [TT] and from a point of view of shift radix system in pQ . 

Conjecture . For all fixed —2 < A < 2, all integer sequences (a n ) defined 
by 

(1) < a n+2 + Aa n+ i + a n < 1 

with initial value (ao,ai) £ Z 2 are periodic. 

For (x,y) = (a n ,a n+ i), we have 

(a n+1 , a n+2 ) = (y, ~[x + Xy\), 

and it defines a map F : (x, y) i— > (y, — [x + AyJ ) on Z 2 . In other words, we 
are interested in the dynamics F on Z 2 : 

G) - (- , -\) © + 

where (x) = x — [x\ . Since the eigenvalues of the matrix are two conju- 
gate complex numbers of modulus one, this dynamics can be regarded as a 
rotation having invariant confocal ellipses, acting on the lattice 1? . After 
the rotation of angle 9 with A = — 2cos#, we translate by a small vector to 
make the image lie in Z 2 . An affine equivalent formulation using Euclidean 
rotation is found in the next section. From the shape of the inequality in ([T]), 
the dynamics ([2]) is reversible, i.e., F(x,y) = (y,z) implies F(z,y) = (y,x). 
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Thus we have c^F^ 1 = F<p with <j>(x,y) = (y,x) and F is a bijection on Z 2 . 
In other words, F is a composition of two involutions F(f> and (j). 

The conjecture is supported by numerical experiments [241 [I]. It is 
also expected from a heuristic ground: cumulation of errors of F n from 
the exact n6 rotation is expected to be small and seemingly impossible 
to avoid hitting the same lattice points. The cumulative error bound is 
discussed in |18|. 126] . However this problem is notorious, and our knowl- 
edge is limited. We only know that Conjecture 1 holds for 11 values A = 
0,±l,(±l±\/5)/2,±\/2,±\/3, see EH IB- Apart from three trivial cases 
0, ±1, the proof is highly non trivial and uses the self-inducing structure 
found in the associated planer piecewise isometry when 6 fix is rational and 
A is quadratic. If 6/tt is rational, then we can embed the problem into piece- 
wise isometry acting on a certain higher dimensional torus (see |20|. [17] , and 
also [U [5] for connection to digital filters) . Piecewise isometries have zero 
entropy |12| . but we know little on their periodic orbits |14j . It is noteworthy 
that a certain piecewise isometry generated by 7-fold rotation in the plane 
is governed by several self-inducing structures [191 021 E] , but it is irrelevant 
to the map F. If A is a rational number whose denominator is the power of 
a prime, then the dynamics is understood as the composition of p-adic rota- 
tion and symbolic shift in [ID], but it seems difficult to extract information 
on periodic orbits through this embedding. At this stage, we are interested 
in giving a non trivial general statement for this dynamics. In this note, we 
will show 

Theorem 1. For all fixed A 6 (—2,2) there are infinitely many periodic 
orbits of the dynamics (0) on Z 2 . 

More precisely, we prove that there are infinitely many symmetric periodic 
orbits (see §3 for the definition). See Corollary [3] in £j6j for a qualitative 
statement. Theorem[T]is new for all A except the above 11 values, and gives 
another support of the conjecture. Note that the idea of observing sym- 
metric periodic orbits dates back to G. Birkhoff, who showed the existence 
of infinitely many symmetric periodic orbits for the restricted three body 
problem (7J [8], whose dynamics is composed of two involutions. 

Adding a counting technique of lattice points in number theory, we can 
generalize Theorem [I] to the sequences generated by: 

—f] < a n+ 2 + Aa n+ i + a n < 1 — 77 

for a fixed r/ 6 R, as in Theorem [2J It covers a significant class of 
discretized rotations, but we do not know how large is this two parame- 
ter family within the set of all invertible discretized rotations in Z 2 , up to 
conjugacy. 

We say p = p(x, y) > is the period of (F n (x, y)) n gz, if it is the smallest 
positive integer p with F p (x,y) = (x,y). If there is no such p, then p(x,y) 
is not defined. It is remarkable that the distribution of periodic orbits dras- 
tically changes by whether 6/tt is irrational or rational. We have 
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Lemma 1. Let 9/tt be irrational and p is a positive integer. Then there are 
only finitely many periodic orbits of period p. 



This fact follows from Theorem 2.1 of [24J but we give a quick proof in §5. 
On the other hand, if 0/tv is rational, in view of the above torus embedding, 
it is natural to obtain infinitely many periodic orbits of period p, which fall 
into the same period cell. Theorem 4.3 in p] gives a concrete example of 
infinite periodic orbits of period p for 9 = (1 — 1/p)tt with an odd prime p. 

By Theorem Q] and Lemma Q3 we know that there exist arbitrary large 
periods, if 9/tt is irrational. We expect the same to hold for all A / 0, ±1, 
but there are proofs only for the above eight quadratic cases. 



2. Setting and strategy 

Let A = —2 cos 9 where 9 is a real number in (0, tt) and let Q 
Our transformation on Z 2 : (x, y) i— > (X, Y) is written as 



sin v cos I 
1 



with pi 6 [0, 1). Since 



3--. ->)(>(: 



cos 9 — sin9\ I J 



Q ' sin0 cos9 ) Q \ - I -A 



we view this algorithm as 



Q' 1 


G 






- esc 6 


) 








cot 9 



Thus it is the dynamics acting on the lattice C = [ n ) Z + ( | 

written as the composition of the Euclidean rotation of angle 9 followed by 
a small translation 

/ cot 6» N 

VH>V + /il j 

with /j, £ [0, 1). Let i? be a positive real number and B(R) be a ball of 
radius R centered at the origin. Define a trap region T(R) by 

'cot 6> N 
1 

The situation is demonstrated in Figure [TJ 

Now we explain the strategy of the proof. It is clear from the description of 
the dynamics, that every unbounded orbit starting from a point in CDB(R) 
must visit at least once the trap region T(R). Assume that there are only 
finitely many periodic orbits of ([3]). Since we are dealing with dynamics on 
the lattice C, periodicity of an orbit is equivalent to its boundedness. We 
argue by contradiction: the assumption that there are just a finite number 
of periodic orbits would imply that, for R 3> 1, the number of lattice points 



T(R) = {x + y 



x e B(R),y e [0,1) \fl(fl). 
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Figure 1. Trap Region 



in T(R) are strictly less than the number of symmetric unbounded orbits 
starting from B(R), which is impossible. 

3. Lower bound of unbounded orbits 

Symmetric periodic orbits of time-reversal dynamics had been studied for 
a long time. We shall make use of a well-known property which may belong 
to a folklore. Let (a n ) be a bi-infinite integer sequence and b be an integer. 
We say that (a n ) is periodic, if a n+b = a n hold for all n. The sequence (a n ) 
is symmetric at b/2, if db- n = a n holds for all n, and it is doubly symmetric 
if it is symmetric at b\/2 and 62/2 with b\ ^ 62- 

Lemma 2. If a sequence is doubly symmetric then it is periodic. Moreover 
symmetric and periodic sequences are doubly symmetric. 

Proof. Assume that (a n ) is symmetric at 61 and 62 with b\ ^ bi- Then 
a„ + 6 2 _6 1 = a bl - n = a n . Let (a n ) be periodic of period b and symmetric at 
c. Then a c _ b _„ = a n+b = a n . □ 

In fact, this lemma is a restatement of Theorems 1 and 2 in [13], where 
symmetric periodic orbits of time reversal dynamics composed of two invo- 
lutions are studied. The reader finds precise description on such periodic 
orbits there. 

Let (x,y) G Z 2 . To the bi-infinite orbit (F n (x,y)) n& % we can associate 
uniquely the bi-infinite sequence (a n ) consisting of the 1-st coordinates of 
the elements of the orbit. It is clear that (F n (x, y)) is periodic if and only if 
(a n ) is periodic. Hereafter we identify the orbit (F n (x, y)) and the bi-infinite 
sequence (a n ) and say that an orbit (F n (x,y)) is symmetric if (a n ) is so. 
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Assume that (a n ) is symmetric at 6/2. If b is odd, then a^_iy 2 = a (6+i)/2 
and the orbit is of the form: 

■ ■ ■ ,C3,C2,Ci,X,X,Ci,C 2 ,C 3 ,. . . 

with some X S Z and a sequence (c n ) C Z. Clearly (c n ) is determined by 
X. Such orbits are in Fix(</>), the set of orbits fixed by the involution eft. If 
b is even, then a&/2-i = a fe/2+i anci the orbit is of the form 

• • • ,C3,C2,Cl,X, Y, X, Cl,C2,C3, . . . 

for some X, Y E Z and (c n ) C Z. Of course (c n ) is determined by X and 
Y. These orbits belong to Fix(F(p), the set fixed by the other involution 
F<p. By Lemma[2j each periodic orbit of (F n (x,y)) is doubly symmetric. It 
belongs to the same fix point set if the period is even, and to the different 
fix point sets if the period is odd. 

Remark 1. Not all orbits are symmetric. For e.g., if A = (1 + y/b)/2 then 
we have 

(-1,4) -> (4,-6) -> (-6,5) -> (5,-3) -> (-3,-1) -> (-1,4). 

We do not know a way to estimate from below the number of asymmetric 
orbits. 

By our assumption of reductio ad absurdum, there exist only a finite 
number, say C±, of periodic orbits. For any R, the number of points in 
£ n B(R) whose orbits are periodic is less than C\. 

3.1. Unbounded orbits in Fix(</>). Let (a n ) be an unbounded orbit in 
Fix(</>). By Lemma[2l (o n ) can not be doubly symmetric, which implies that 
two orbits starting from different fixed points in C never intersect. Apart 
from a finite number of exceptions, points of the form 

/— csc#\ „ (cot 9\ ^^t~,/t-,\ 

V ) v 1 ) €CnB ^ 

generate distinct unbounded orbits. From 

X 2 (- esc 9 + cot O f + X 2 < R 2 , 

we conclude that there are at least 2i?cos(#/2) — C\ unbounded orbits in 
Fix(</>) starting from C f] B(R). 

3.2. Unbounded orbits in Fix(i ? 0). Similarly, by LemmaO (a n ) can not 
be doubly symmetric. So, by (pQ), our task is to count the number of the 
pairs (X, Y) which satisfy 

(4) < X + XY + X < 1 
and 

f— csc8\ TA /cot^\ „ /T , N 

(5) Xl Q 1 ) <=B(R). 
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For this computation, we substitute the inequality (JH) by 

(6) - 1 < X + XY + X < 1 

and count the number of pairs (X, Y) satisfying ([6]) and It is clear that 
for a fixed Y, there is a unique X which satisfies ([6]). Since X = Y cos 8 + e 
with |e| < 1/2, we have 

(7) (Ycot#-Xcsc0) 2 + Y 2 = Y 2 + — ^ < R 2 

sin 

from ([5]), and we have at least 2R — C2 such points with a non negative 
constant C2. 

If A is irrational, then there is no (X, Y) 7^ (0, 0) which satisfies either 

-l = 2X + XY or = 2X + AY. 

Using the symmetry (X, Y) <-> (—X, — Y"), we see that the number of (X, Y) 
with ([5]) and 

-1<X + XY + X<0 

is exactly one less than the number of (X, Y) with dH) and ©, which counts 
the origin. Thus the number of (X, Y) having (jj]) and ([5]) is at least R—C2/2. 

If A is rational, then we additionally have to take care of the points (X, Y) 
on the line — 1 = 2X + XY and = 2X + XY. However we can easily show 
that the number of (X,Y) with (J5J on the line -1 = 2X + XY and the 
one on the line = 2X + XY differ only by some constant. Thus in any 
case, there are at least R — C3 unbounded orbits in Fix(F(p) starting from 
C n B(R) with a non negative constant C3. 

4. Lattice points in the trap region 

By construction, if the trap region T(R) and the line x = ycot6 + c has 
non empty intersection, then it is a half-open interval of length csc#. Thus 
if the line x = y cot 6 + c intersects C n T(R) then it is a single point. The 
lattice C is covered by a family of parallel lines: 

E = {x = y cot 8 — k esc 8 \ k £ Z}. 

We easily see that the distance between adjacent lines of S is 1. Thus there 
are exactly 2[R\ + 1 points in £ n T(R). 

Remark 2. If 6 < 2vr/3, then 2Rcos{8/2) - d + R - C 3 > 2R+ 1 holds for 
sufficiently large R and we immediately obtain the desired contradiction. 

Let us take into account the reversibility of F. Since we are dealing with 
unbounded symmetric orbits starting from £ n B(R), if (o n , a n +i) = (C, D) 
then there is an index m such that (a m , a m+ \) = (D, C). Let <3? : £ 1— > £ be 
defined as follows 

f-csc8\ fcot9\ f-csc8\ /cot0\ 

J+4 1 J^H ° ) +x \ 1 J" 
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If an orbit visits <3?(T(it!)) n T(R) then the number of visits is at least two. 
In other words, we only have to count the number of lattice points up to 
this symmetry by $ in T(R). 

— csc ~\~ cot 

because the two vectors f C ^ C ^ and ^ C °^ ^ have the same length. Thus 

the reflection leaves the vector f C os(0/2)^ = f ^ ^ 1 ^ ^1 cos(9/2) 
invariant. 

To make computation easy, we rotate T(R) and C by —6 and present the 
situation in Figure [2j T(R)', £ are the images by this rotation and \& is 



¥(T(R» 




Figure 2. Symmetry of the trap region 

the corresponding reflection. Then every line y = k with k G 1*P\[—R,R] 
contains a single point in CJ n T(R)' and the reflection \E' leaves the vector 

/ cos# sin#\ /-sin(0/2)\ _ /sin(0/2)\ 
\^-sin0 cos0y V cos(0/2) J ~ \cos(9/2)J 

invariant. 

So our task is to estimate from above the number of lattice points in 
C n T(R)' which are below the line y = a: cot (0/2). The intersection of the 
line and the boundary of T{R)' with the largest y-coordinate is 
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and we have 

«* G) - 1 + 5 = (§) + 5 + ° (s) ■ 

We count the number of points whose y-coordinate does not exceed this 
value, i.e., the points in the shaded part in Figure [2j Thus the number 
of lattice points up to symmetry in £ n T(R) is bounded from above by 
R + Rcos(9/2) + C4 with a non negative constant C4. 



5. Proof of Theorem CD and Lemma Q] 

From the assumption that there are only finitely many periodic orbits, we 
derived several estimates in the previous sections. By Lemma El unbounded 
orbits in Fix(0) and those in Fix(FcJ)) have no intersection. Thus 

2Rcos{9/2) -d + R-C 3 

distinct unbounded orbits must visit T{R) and there are only 

R + Rcos{9/2) + C 4 

lattice points in C fl T(R) up to symmetry. However 

2Rcos(9/2) -d + R-C 3 <R + Rcos{9/2) + C 4 

does not hold for sufficiently large R. The proof of Theorem [T] is finished. 

Let us show Lemma [TJ Assume that there are infinitely many (x,y) that 
F p (x,y) = (x,y) with p > 2. By induction using (J3|), we have 



where 
if 



COS} 

sin? 



— sm'j 
cos p> 



6 C and ||vj|| < esc 9. Here 
is sufficiently large, then 



cos p9 
sin f)9 



- sm-j 
cos p 



Here c is computed as the operator norm: 

1 



i=l 



is the Euclidean norm. However 



> c 



1 



- cos p9 
sin p9 



sin p9 
1 — cos z 



-1 



sm 



> pcsc9. 



p9_ 
~2 



which is positive since 9/ir is irrational. This gives a contradiction. 
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6. Generalization 

One can generalize the result to the sequences defined by: 

-rj < a n+2 + \a n +i + a n < 1 — rj, 

with rj € R. This kind of interval shifts are studied, for e.g., in |17[ 12 1 J . In 
complete analogy to our main result, we have 

Theorem 2. For a fixed A G (—2,2) and r/ £ R, there are infinitely many 
periodic orbits of the dynamics 

(8) {x,y) ^ {y,-[Xy + x + r]\) 
on 1?. 

We can deduce a qualitative statement: 

Corollary 3. There is a positive constant C depending on A and 7] such 
that within B(R), there are at least CR periodic orbits of |5)/ 

Hereafter we sketch the proof of Theorem [2] and Corollary [3j Putting 
n = 77/(2 + A), the inequality becomes 

< (a n+2 + k) + A(a n+ i + k) + (a n + k) < 1. 

Therefore by substituting C with £' = C + Q^ 1 , our algorithm has 

exactly the same shape as (|3|). Though the error term becomes worse than 
the one in §4, we can show that the number of lattice points of £ up to 
symmetry within the trap region is 

R + Rcos(9/2) + 0(R 2/3+e ) 

for any positive constant e. Here we used the method of Vinogradov to count 
the number of lattice points in the cylindrical region bounded by curves of 
positive curvature, for e.g., see p. 8-22 of [23] or [22| [16] . 

Similarly to §3, there are 2Rcos(9/2) — C\ unbounded orbits in Fix(</>). 
We count unbounded orbits in Fix(i ? ^), i.e., the number of (X, Y), which 
satisfy: 

(9) AF/2 mod in [-77/2, (l-r/)/2) ^0 
and ©. 

If A is irrational, then (A/2)Y mod 1 is uniformly distributed and the 
number of such Y's is R + o(R). When A is rational, put A/2 = p/q with 
(p,q) = 1. Then (A/2)Y = i/q mod 1 for i £ {0, 1, ... ,q— 1} with the same 
frequency 1/q. 

Let us study the case that q is even. Since 

{i/q mod 1} n [-??/2, (1 - r/)/2) 

has cardinality q/2, the number of points with Q and ([7|) is again R + o(R). 
Once we have this estimate R + o(R) then 

2i?cos(#/2) - Ci + R + o{R) < R + Rcos(6/2) + o(R 2/3+e ) 
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does not hold for sufficiently large R. Here C\ appears only once in the left 
side, because it is the number of periodic orbits of the system. We obtain 
the contradiction for R 3> 1. 

It remains to show the case when q is odd. Then 

{i/q mod 1} n [-r?/2, (1 - rj)/2) 

has cardinality either (g— 1)/2 or (q+l)/2 depending on r/. Thus the number 
of unbounded orbits in Fix(Fcj)) is bounded from below by R — R/q + o(R). 
Thus we have to show that 

2Rcos(9/2) - d + R - R/q + o(R) > R + Rcos{9/2) + o(R 2 / 3+e ) 

for large R. This is valid because 

cos(0/2) > l/q 

holds for q > 2, since cos(0/2) = ^(1 + cos(0))/2 = y^l -p/q)/2. 

The validity of statement of Corollary [3] is invariant under affine trans- 
formations fixing the origin, up to appropriate changes the constant C. To 
prove Corollary [3l note that we may take 

Ci = (1 - e)22 (cos(0/2) - l/q) 

with a small e > in the above proof to get the same contradiction. Here 
the term —l/q is necessary only in the last case of the proof. 
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